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Abstract 
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CSJ ■ We previously proposed a formula for calculating the Coulomb energy of spherical nucleus with 

jrt ■ Wood-Saxon charge distribution. In this work, the analytical formula is extended for description 
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of the Coulomb energy of nucleus with /32 deformation. 



Corresponding author : lium_816@hotniail.com 



The calculation of the Coulomb energy for complicated charged system with small com- 
puting effort and high accuracy is a great challenge in physics and quantum chemistry 
-3|]. For a system with an arbitrary charge distribution p(r), the direct term of 



research 



the Coulomb energy can be calculated with 

2 J J |r-r'| 

However, the six- dimensional integration in Eq.(l) is very time-consuming and becomes a 
bottleneck in the large-scale calculations of potential energy surfaces of nuclear systems. In 
this work, we attempt to propose an analytic expression for calculating the Coulomb energy 
of nucleus with both the nuclear surface diffuseness and nuclear P2 deformation being taken 
into account. For reader's convenience, the approach to calculate the Coulomb energy of a 
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nuclear system proposed in our previosuly work [l[ is reviewed firstly, and then the analytical 
formula for calculating the Coulomb energy of spherical nucleus with Woods-Saxon density 
distribution will be extended for description of the Coulomb energy of nucleus with (32 
deformation. 

The Coulomb energy of an arbitrary nuclear system can be obtained by 

Ec = lj p{v)Vc{v)dv, (2) 

where Vc'(r) is the Coulomb potential which is obtained by solving the Poisson equation 

VVc(r) = -47rep(r). (3) 

The charge distribution of a nucleus is usually described by a Woods-Saxon form, 

l + exp(^) 

Where, po and a denote the central charge density and the surface diffuseness, respectively. 
TZ defines the distance from the origin of the coordinate system to the point on the nuclear 
surface. For an axially deformed system, TZ is expressed as, 

n{e) = R^[i + hY2o{e) + -]- (5) 

In the calculation of Coulomb energy of nucleus as a function of nuclear deformation, we 
remain the central charge density po of the nucleus unchanged by using the conservation of 



charge number and varying the half-density radius i?o to consider the effect of incompress 
ibihty of nuclear matter in the nucleus. 



fl 



The Poisson equation is solved by a code hwscyl (a Fortran subroutine in FISHPACK j^ ) 
which is an adaptive fast solver for solving a five-point finite difference approximation to the 
modified Ifelmholtz equation in cylindrical coordinates using a centered finite difference grid. 
We calculate the Coulomb potential in cylindrical coordinates within a region x = ~ 40 
fm and z = —40 ~ 40 fm (using a grid with step size 0.1 fm). It is known that when r ^ TZ, 
the asymptotic behavior of the Coulomb potential of a nucleus is Vc = eZ/r, which gives the 
boundary condition in solving the Poisson equation. The Coulomb energy of an arbitrary 
axially deformed nuclear system can be obtained with a two-dimensional integration over the 
Coulomb potential Vc{r) which can be calculated with the very fast solver for the Poisson 
equation mentioned above. 

In our previous work [l|, we investigated the Coulomb energies of spherical nuclei with 
Wood-Saxon charge distributions. The central charge density of a nucleus is obtained with 
the Skyrme energy density functional together with the extended Thomas-Fermi (ETF) 
approach [5|. The nuclear surface diffuseness a varies from 0.1 to 1.2 fm in which the central 
charge density is remained unchanged. We find that the Coulomb energies of spherical nuclei 
with Woods-Saxon charge distributions can be well described with an analytical expression 
based on the leptodermous expansion |6[, 

Ecoui = EfF{uj) (6) 

with (l| 

5 

F{uj) = 1 - -w^ + c^uo^ + uj^ + c^u^ + CqU^ + ... (7) 

Where, '^ = -%% and E)-, denotes the Coulomb energy of a spherical nucleus with uniform 
charge distribution, 

4"' - ^^. (B) 

Z denotes the charge number of the nucleus and R = [Z/(^po)] is the corresponding 
radius of a spherical nucleus with uniform charge distribution. By fitting the calculated 
Coulomb energies with numerical integrations for a number of spherical nuclei along the 
/3-stability line, we obtained the coefficients C3 = 3.005, C5 = —4.822, cq = 2.934. 
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FIG. 1: (Color online) Coulomb energies of ^^^U as a function of nuclear surface diffuseness a and 
quadrupole deformation /32 (a); and as a function of /32 and fS^ (b). 

With the same approach, we investigate the Coulomb energies of deformed nuclei. In 
Fig.l (a), we show the calculated Coulomb energy of ^ss-q ^^ g^ function of nuclear surface 
diffuseness a and quadrupole deformation P2, and we show the corresponding Coulomb en- 
ergy as a function of (32 and (3^ (with a = 0.55 fm) in Fig.l(b). One can see that the Coulomb 
energy decreases with increase of the nuclear surface diffuseness and of the deformation. In 
this work, we write the Coulomb energy of a nucleus as 



Ecoui = E^^^F{co)G{u,(3), 



(9) 



with a factor G{uj, (3) to consider the influence of nuclear deformation. For nucleus with (32 
deformation, we assume that the factor G{uj, (3) has a form 

1 

471 '■ 



G{U,P2) 



-(3l + hiuj(3l + h2U^(3l + h(3l + 64/32 + - 



(10) 
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FIG. 2: (Color online) (a) Coulomb energy of 238]j g^g g^ function of (^2 deformation (with a = 0.55 
fm). The squares denote the results with numerical integration [with Eq.(2)]. The dot-dashed 
curve and the solid curve denote the results of Eq.(9) without and with the 6i to 64 terms being 
taken into account, respectively, (b) Relative deviations |£'coui — Ec\/Ec of the Coulomb energies 
from the numerical integration results, ii'coui denotes the calculated Coulomb energy with Eq.(9) 
and (10). The region with light gray denote that the relative deviations are smaller than 0.05%. 



TABLE I: Values of w = -^-^ for some nuclei by taking a = 0.55 fm and R = 1.2j4^" fm. 



16 



O 



40 



Ca 



90 



Zr 



144 



Sm 



208 



Pb 



238 



U 



298 



114 



CJ 



0.33 0.24 0.19 



0.16 



0.14 



0.13 



0.12 



Where the term —-^(^l is presented by Greiner and Maruhn in Ref. [TJ. By fitting the 
calculated Coulomb energies with Eq.(2) as a function of surface diffuseness (a < 0.7 fm) 
and quadrupole deformation (|/32| < 0.5) for a number of nuclei along the /3-stability line, 
we obtain the coefficients hi = ^, &2 = 0.188, 63 = —0.007 and 64 = 0.018. 

In Fig.2(a), we show the Coulomb energy of ^ssij ^^ g^ function of nuclear quadrupole 
deformation /32. The squares denote the results with Eq.(2). The dot-dashed curve and the 
solid curve denote the results of Eq.(9) without and with the hi to 64 terms being taken 
into account, respectively. One can see that the higher-order terms of deformation are still 



required for system with large deformation. In Fig.2(b), we show the relative deviations 
l-E'coui ~ Ec\/Ec of the Coulomb energies from the numerical results for a number of nuclei 
A = 16 ~ 300 varying the surface diffuseness (a < 0.7 fm) and the quadrupole deformation 
(I/32I < 0.5). -Ecoui denotes the calculated Coulomb energy with Eq.(9) and (10). From 
Fig.2, one can see that the Coulomb energy obtained with the analytical formula Eq.(9) is 
close to the calculated results with numerical integration for most cases. In Table 1, we 
list some typical u values of a series of nuclei from light to heavy. For intermediate and 
heavy nuclei, u has a value about 0.1 ~ 0.25, the corresponding relative deviations of the 
Coulomb energies with Eq.(9) are smaller than 0.05% [denoted by light gray in Fig. 2(b)] for 
almost all cases with a < 0.7 fm and |/32| < 0.5. For light nuclei, the corresponding values 
of u are larger than 0.25 in general and the relative deviations of the Coulomb energies 
with Eq.(9) slightly increase for some cases with strong deformations. It is known that the 
charge distributions of light nuclei are usually described by gaussian functions rather than 
the Woods-Saxon form. The analytical expression of the Coulomb energy of a system with 
gaussian charge distribution can be found in Refs. [ll, |6|. 

In summary, the Coulomb energy of axially deformed nucleus with Wood-Saxon charge 
distribution has been investigated. The Coulomb energy of a nuclear system was numeri- 
cally calculated with a two-dimensional integration over the Coulomb potential which was 
obtained by solving the Poisson equation. By fitting the numerically calculated Coulomb 
energies for a number of nuclei from A = 16 to 300 with P2 deformed Wood-Saxon charge 
distribution, an analytical formula, that is a function of nuclear /32 deformation and surface 
diffuseness, is finally obtained. The relative deviation of the Coulomb energy with the pro- 
posed formula is generally smaller than 0.05% for intermediate and heavy nuclei at normal 
deformations. For nuclear system with |/32| > 0.5, such as fissioning system, the proposed 
formula could not be applicable, and the two-dimensional numerical integration over the 
Coulomb potential has to be performed to obtain accurate results. 
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